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Abstract. Networks of random quantum scatterers (S-matrices) form paradigmatic models 
for the propagation of coherent waves in random media. S-matrix network models cover uni- 
versal localization-delocalization properties and have some advantages over more traditional 
Hamiltonian models. In particular, a straightforward implementation of real space renormal- 
ization techniques is possible. Starting from a finite elementary cell of the S-matrix network, 
hierarchical network models can be constructed by recursion. The localization-delocalization 
properties are contained in the flow of the forward scattering strength ('conductance') un- 
der increasing system size. With the aid of 'small scale' numerics qualitative aspects of the 
localization-delocalization properties of S-matrix network models can be worked out. 

Keywords: S-matrix network models, Localization-delocalization transition, renormaliza- 
tion group, multifractality 



1 Introduction 

Localization and fluctuation phenomena are fascinating universal features in any dis- 
ordered coherent wave mechanical system. In 1982 Shapiro Q pointed out that a 
convenient modeling of disordered coherent wave mechanical systems can be given 
by networks of unitary random scattering matrices. This modeling covers essential 
symmetries and characteristic length scales, but docs not rely on particular disper- 
sion relations and specific details. Therefore, it is of relevance in various fields of 
theoretical physics, e.g. in optics, mesoscopic electronics, and quantum chaos. S- 
matrix network models (NWMs) have a number of advantages over more traditional 
Hamiltonian models. NWMs yield directly transport quantities ||, wave packet dy- 
namics ||, and quasi-energy eigenvalues and eigenstates A well known NWM 
introduced by Chalker and Coddington (CC-model) Q (see also describes the 
situation of disordered two-dimensional (2D) electrons undergoing the quantum Hall 
localization-delocalization (LD) transition in a strong perpendicular magnetic field. 
A major advantage of S-matrix NWMs is a straight forward implementation of real 
space renormalization group (RG) techniques [| which open new perspectives for 
investigating LD transitions. 

Here we address the construction of S-matrix NWMs starting from symmetry consid- 
erations and/or from an Anderson tight-binding Hamiltonian. We discuss alternative 
implementations for a RG action. To reach a qualitative overview of phases in the 
phase space of scattering strengths one can construct hierarchical network models for 
which a RG algorithm can be formulated that is designed for 'small scale' numer- 
ics. The construction of S-matrix NWMs and the RG methods are demonstrated on 
particular two-dimensional models. 



2 Construction of S-matrix networks 

Quite generally, a NWM can be constructed as follows. Take a regular network of 
N sites and N bonds. Each bond a carries propagating wave modes (n+ incoming- 
modes and n~ outgoing modes) represented by complex amplitudes, ip n ±- On the 
sites unitary S-matrices map incoming to outgoing amplitudes. The elements of each 
S-matrix are (in general) random quantities respecting the symmetries of S and are 
characterized by typical scattering strengths. Random phases are attached to the 
amplitudes on the bonds. They simulate the random distances between scatterers in 
realistic systems. 

2.1 Topology and site S-matrix 

The construction of a NWM is fixed by the choice of a certain type of random S- 
matrix and a J\f x TV connectivity matrix C which has elements Cjj = 1 if a wave 
mode can propagate from site i to site j and CL- = otherwise. The NWM defines a 
dynamical system by S and C. Together they lead to a unitary matrix U [|[ ^, (To) 
called 'network operator'. It maps all incoming to outgoing wave modes. We choose 
a convenient time unit and denote the vector of the B cx N bond amplitudes as ip 
such that ip(t + 1) = Utp(t). The eigenphases cj> n (n = 1 . . .B) and corresponding 
eigenvectors ipn of U, Uip n = exp(i</j n )-0 n , can then be interpreted as quasi-energies 
and corresponding eigenstates. 

2.2 Two-dimensional networks 

Let us now report on two-dimensional NWMs with the topology of a square lattice 
and a connectivity matrix that connects each site with its four nearest neighbors. The 
corresponding S-matrix is graphically represented in Fig. 1 and reads 

/ r 1 4 1 t 31 
df r 2 4 2 
* - t 13 df r 3 

The parameters of each S-matrix are the amplitudes of transmission t, reflection r, 
and deflection to the left dh and right dn which contain random phases compatible 
with symmetry requirements. Scattering strengths of a site S-matrix allow to define 
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Fig. 1: U2 network model with unitary S-matrices. On the left the elementary S-matrix 
and on the right a section of the whole network are graphically represented. 
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transport mean free paths l e (measured in units of the lattice constant) for a prescribed 
forward direction by l/l e = (1 —T)/T. Here T is the probability to be transmitted 
in the prescribed forward direction. Such definition is known from the scattering 
theoretical approach to quasi-lD systems (e.g. p"T[). 

The unitarity of the S-matrix requires | £ | 2 + 1 r | 2 + 1 | 2 + 1 | 2 = 1. The corresponding 
parameter space of scattering probabilities, a quarter of a pyramid, is shown in Fig. 2 
and forms the phase space of the U2-model ('IP stands for unitary, '2' for the two- 
dimensional lattice structure). The points in the corners of the phase space are the 
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Fig. 2: The three-dimensional phase space of the U2- model is spanned by scattering 
strengthsO < \t\ 2 , |r| 2 , \d L \ 2 < 1 with \t\ 2 + \r\ 2 + \d L \ 2 = l-\d R \ 2 and \d R \ 2 < 1. The shaded 
area corresponds to the U2NC-model for which \d L \ = \d R \. The points LHFP, RHFP, LFP 
and MFP are characterized by = 1, \d R \ — 1, \r\ = 1 and \t\ = 1, respectively. They 
form the corners of the U2-phase space. The critical fixed point of the Chalker-Coddington 
model (CCFP) corresponds to \d L \ 2 = \d R \ 2 = 1/2. 

metallic fixed point (MFP) (\t\ = 1), the localization fixed point (LFP) (|r| = 1), and 
the left (right) handed fixed point LHFP (RHFP), \d L \ 2 = 1 {\d R \ 2 = 1). The term 
'fixed point' refers to the renormalization group action to be discussed below. 

2.3 Symmetries and phase space 

In the absence of handedness of the scattering process, i.e. \di\ = \dn\, the model 
is denoted as non-chiral (U2NC-model). The phase space of the U2NC- model corre- 
sponds to the shaded plane in Fig. 2. The point of \d^\ 2 — \dn\ 2 = 1/2 is denoted 
as Chalker-Coddington fixed point (CCFP) since the U2-model reduces to two un- 
coupled CC-models for vanishing t and r and the phase space of the CC-model is 
one-dimensional (ID) and parameterized by the values of T ;= |g£l| 2 (see Fig. 3a). 
Further symmetries of the S-matrix lead to a reduction of the number of independent 
parameters. This typically eliminates random phases, but can also lead to restrictions 
in the phase space of scattering strengths. For time reversal symmetric scattering in 
the U2NC-model the S-matrix is symmetric and the phase space is reduced by \r\ + 
\t\ > 1. The model is then denoted as 02NC-model where 'O' stands for orthogonal 
(the network operator can be diagonalized by orthogonal matrices). Its phase space 
is shown in Fig. 3b. Imposing symmetry constraints reflecting a sublattice structure, 
e.g. SaS'Sg = S\ where E 3 := diag(l, 1, -1, -1), the U2NC-model becomes 'S 3 - 
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Fig. 3: a) The phase space of the Chalker-Coddington model. Two of the fixed points cor- 
respond to quantum Hall plateaus (LHFP and RHFP), the critical point CCFP corresponds 
to the plateau to plateau transition, b) The phase space of the 02NC-model in \r\, \t\ pa- 
rameter space. The time reversal symmetry leads to the restriction \t\ + \r\ > 1 (grey area), 
c) The phase space of a chiral model. The chiral symmetry (see text) leads to the restriction 
|i| 2 > k| 2 (grey area). 



chiral'. Although \dz\ — \dn\ is still valid, the chiral symmetry leads to further phase 
space constraints, e.g. to \r\ 2 < \t\ 2 (for d ^ 0) for E3-chirality (Fig. 3c). 

2.4 From tight binding models to S-matrix networks 

So far S-matrix network models appeared as models on their own, dictated by topology 
and symmetry. However they can also appear as the result of a mapping from a 
microscopic Hamiltonian, as is the case for the CC-model ||. In this subsection we 
show how tight binding models can also lead to S-matrix network models. This opens 
a new flexibility in the methods to study these models. 

The general idea is to take a typical microscopic cell of a given tight binding model 
and attach ID semi-infinite leads to it, as many as is appropriate for topology and 
symmetry aspects. Then calculate the S-matrix of this device. It may then serve to 
define a NWM. 

We proceed as follows. The Hamiltonian H = Hq + Hi is decomposed into a part 
that describes only the leads, Hq, and the part that describes the microscopic cell 
together with its coupling to the leads, H\. The scattering states with continuous 
spectrum are chosen to obey energy normalization (a(E)\a(E')) = 5(E — £"). Then 
the S-matrix is given as S a p — 5 a p — AT{iT a p where the T-matrix (operator) is related 
to the resolvent (Green's function) G+ = (E + i0+ - H)- 1 by T = H x + H X G + H\ . 
Note that Hq describes 'isolated' terminals, i.e. (semi-infinite leads) and the modes 
| a) contain both 'in'- and 'out'-states. Introducing the projector Q = 1 — P onto the 
Hilbert-space of the terminal states the T-matrix elements read 

T ali = ((3 | HqpGppHpQ | a) . (2) 

The projected Green's function Gp P can be determined by a non-hermitian effective 
Hamiltonian defined entirely on the P-space, i.e. we have to invert a V dimensional 
matrix when V is the number of sites in the microscopic cell. 
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It is convenient to treat the terminals as semi-infinite ID tight-binding models with 
only kinetic energy and unit hopping strength. The terminal Green's function as well 
as the amplitudes of the terminal states can be calculated analytically. Together with 
the solution of the V x V matrix problem the S-matrix can be obtained from Eq. (0) . 
To make the discussion more explicit we present the results for two models Of. 
First consider the microscopic cell in Fig. 4a consisting of a on-site random energy e 
attached to four terminals. For energies in the band center of the ID terminals the 
scattering matrix elements are |b| 

8 - 2ie 

t = d L =d R = — - — 2 
16 + e z 

and r = 1 — t. The model is a member of the 02NC-model since S is symmetric. 
To get scattering strengths characteristic for an ensemble of disorder realizations one 
should average them with respect to an appropriate distribution of e. 
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Fig. 4: Elementary S-matrix for tight binding Hamiltonian models: a) A central site with 
random on site energy e is coupled to four semi-infinite leads (terminals) characterized by 
unit hopping amplitudes, b) Four central sites enclosing a random flux <j> are coupled to 
semi-infinite leads. 



Second we address a microscopic cell for the random flux problem as shown in Fig. 4b. 
No on-site disorder is introduced, but a random flux phase <j> (the flux quantum is set 
to 2tt) is picked up when moving around the four inner sites (anti-clockwise). At the 
band center the resulting scattering amplitudes (for an incoming wave from the left 
in Fig. 4b) are found to be 

r = (4sm 2 (0/2)-l)/D, t = -4cos(cj>/2)/D , 
4 2 (0) = e -^ /4 (2ze^ - 4i)/D = -(4 1 )* 

with D = 4sin 2 (0/2) + 5. A close inspection shows that the model is 'cr3-chiral', 
i.e. a^Sas — with 03 := diag(l, —1, 1, —1). This property is lost for energies not 
exactly at the band center. 



3 Implementation of a renormalization group 

To detect localized, critical and delocalized phases of an infinite NWM one likes to 
implement a convenient renormalization group (RG). Following the scaling theory of 
Abrahams et al. |l4| one could think of computing the two-probe conductance g{L) 
for a network of size L 2 (Fig. 5a). This involves strong mesoscopic fluctuations in g(L). 
For that reason the self-averaging object A(M) — £{M)/M has been introduced [|l5) 
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Fig. 5: Different ways of implementing a renormalization group: a) The two probe conduc- 
tance g(L) is studied for increasing system size L. b) The finite size scaling variable A(M) 
is studied as a function of width M in quasi- ID geometries, c) The elementary S-matrix can 
be generalized to an S-matrix S(L) for arbitrary system size L. 

to study the scaling behavior. Here £(M) is the localization length corresponding to 
quasi-lD geometries (Fig. 5b) of length L — > oo and finite width M, This so called 
'finite size scaling' method is the 'first choice' method for numerical investigations of 
scaling properties Jul]. It has been successfully applied to the CC-modcl [0, |l7j and 
to the U2NC-model (and variants of it) |jl8|| . Still, it involves 'large scale' numerics, 
i.e. powerful computation units and/or large storage capacities. One reason for this 
is that the finite size scaling is an 'exact' numerical method - no approximations in 
the modeling are involved. 

Alternatively, one can treat the whole system of linear size L 2 as a big scattering 
unit (see Fig. 5c) giving rise to a RG flow of the S-matrix: S a p(l) — ► S a p(L). 
Consequently, this RG procedure allows to study a RG flow in the phase space of 
scattering strengths. In particular the fixed points of this RG flow are of interest. 
The repulsive or attractive behavior of the RG flow close to the fixed points gives 
valuable information about phase transitions. For example, the CCFP is repulsive 
within the phase space of the CC-model while LHFP and RHFP are attractive fixed 
points (see Fig 3a). These simple facts contain essentials of the quantum Hall effect. 

4 Hierarchical S-matrix networks 

The numerical effort in computing a single S-matrix element of S(L) is of the same 
order as in computing a single eigenstate of the network operator U. Furthermore, 
one needs computations for different disorder realizations to take large fluctuations 
into account. Thus, a further idea is needed to make advantage of the RG flow in 
phase space. We give up the 'exactness' of the RG method and consider simplified 
NWMs by recursive constructions of 'hierarchical' NWMs as shown in Fig. 6 for the 
U2-model. A section of a given NWM containing b 2 sites and V 'active' elementary 
S-matrices (S-matrices with only trivial elements 0, 1 are not counted as active) is 
taken. The scatterers are connected such that the resulting network has the same 
terminal structure as each of the elementary S-matrices. Accordingly, the hierarchical 
models are denoted as [b/V]-models. The hierarchical models have fractal dimension 
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Fig. 6: An example (denoted as U2[2/4]) for a recursive construction of a hierarchical 
network model. Four elementary S-matrices are connected such that the resulting 4x4 
network has the same terminal structure as each of the elementary S-matrices. 



= in y/ In 6. In the V -> oo limit the original NWM may be recovered exactly 
(depending on the construction). Obviously, the RG flow works by recursion and 
the computational effort simplifies drastically. The hierarchical NWMs approximate 
the original ones in an uncontrolled manner. Better approximations can only be 
gained by increasing the size of the elementary cell and by studying the convergence 
of quantitative results, e.g. critical exponents. However, qualitatively correct results 
for the existence or non-existence of phase transitions can be obtained already for 
moderate sizes of the elementary cells. The RG method for hierarchical NWMs yields 
qualitative aspects of the RG flow in phase space of the original NWM relying on 
small scale numerics. 



5 Renormalization algorithm 

In the following we outline the six-step algorithm of the RG in hierarchical network 
models. 

(i) In the first step S a p(L = b) is determined as a function of the scattering ma- 
trix elements of all active scatterers for a given elementary cell (respecting unitarity 
constraints). For small numbers V this can be done by hand or computer algebra. 

(ii) The set of scattering strengths T is drawn from an initial distribution,P (T). 
Phases are taken uncorrelated and uniformly distributed in the interval [0, 2ir]. For 
each realization T(L = b) is calculated with the help of the algebraic expression 
obtained in the first step. 

(iii) After collecting a large number (~ 10 3 ) of scattering strengths an approximate 
distribution function Pt=b(T) is determined which depends on the initial distribution 
-Po(T). In particular, it yields the average value (T) (L = b) and other characteristic 
parameters of the distribution. 

(iv) Steps (ii) and (iii) are repeated, but now T is drawn from the distribution Pi = b(T) 
obtained before. This leads to a reasonable approximation for P £=b 2(T) and for 
(T)(i = 6 2 ). 

(v) Step (iv) is repeated N times {N ~ 10 1 ) which yields, for each initial distribu- 
tion -Po(T), the flow of the distribution function Pl(T) and of the mean (T) (L) as 
functions of In L = (N + 2) In b. 
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(vi) Finally, the RG process ((i) - (v)) will be repeated for a number of different initial 
distributions Po(T) to study the RG flow of typical values close to the fixed points. 

6 Small scale numerics 

RG calculations following this algorithm have been performed for the CC-model |6) 
(see also f§) and the Manhattan model (MM) considered in Q. The MM is defined 
by r — and an alternating structure of d^ = or dji = (see Fig. 8) leaving the 
forward transmission T := |i| 2 as the single parameter for its phase space. 
The results for a number of CC[b/V]-models |6) can be summarized as follows: A LD 
transition occurs, characterized by a repulsive fixed point distribution P*{T) (here 
T := |<iL| 2 ) which is broad over the interval [0,1]. After a few iteration steps the 
flow becomes independent of the initial distribution and can be characterized by the 
average value (T). The critical distributions were found to be symmetric around the 
critical mean value T* = 0.5 and developed a shallow minimum there. The critical 
distribution for the CC[3/5]-model and the scaling of its mean value is shown in 
Fig. 7. Furthermore, one observes that with increasing elementary cells the critical 
exponent of the localization length approaches the value of the original CC-network 
in a monotonic way from above. 

In addition to the localization length exponent a second critical exponent is required 
to characterize the LD transition, the critical exponent, a$ — D, of the order parameter 
(the typical density of states) [g(J. This quantity scales to zero at the critical point, 

Pt(L) oc L~( a °- D ) , a Q > D, 

since critical states are not totally space-filling, but multifractal. Here D is the 
dimension of the system which may be fractal in hierarchical NWMs. Adopting the 
method of Klesse and Metzler || to calculate wave functions as eigenstates of the 
network operator U the multifractal exponent ao can be computed for a given network 
model. To reduce the computational effort, all scattering strengths can be fixed to 
the critical value T*. Moderate sizes (L « 50) are used and several realizations to 
obtain ao- It turns out that ao — D^ h / y ^ monotonically approaches the value of the 
CC-model |§ from above as Z)[ b / V ] approaches d — 2 from below. 




# iterations 



Fig. 7: The critical distribution of scattering strengths (left) and the corresponding flow of 
average transmission strengths T under increasing number of renormalization steps (right) 
for the CC[3/5]-model. 
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Fig. 8: The [3/5]- model corresponding to the Manhattan model (top) and the flow of 
the distribution of the forward scattering strength (bottom) under renormalization (L = 
3, 9, 27, 81) to strong localization. Amplitudes on open ends with the same symbol (bullet 
or shaded circle) are identified within one RG step. 

It is believed that a crucial ingredient of the CC-model is its handedness reminiscent of 
the strong magnetic field in quantum Hall systems. To check whether the handedness 
is essential for the occurrence of the LD transition we performed a RG process for 
the Manhattan model which has a one-dimensional phase space like the CC-model, 
but no handedness (even for strong deflection there are as many right circulating 
loops as left circulating loops in a large network) . Indeed, we find no indication for a 
LD transition in the MM. Instead, any initial distribution of the forward scattering 
strength T flows to the strong localization fixed point LFP (see Fig. 8). 

7 Conclusions 

S-matrix network models provide a flexible tool to study disordered coherent wave 
mechanical systems. The construction can be based on topology and symmetry re- 
quirements, but also on microscopic models such as tight binding Hamiltonians. For 
a given network model a recursive construction of associated hierarchical network 
models can be performed. Iterating this construction yields an approximation of the 
original network model which allows to study the flow of scattering strengths in the 
phase space of the network model by a six step renormalization group algorithm. The 
precision of critical parameters obtained by this procedure can be improved by in- 
creasing the number of realizations in each iteration step. The critical parameters are 
then more precise for the particular hierarchical model, however have uncontrolled 
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deviations from that of the original model. Better approximations can only be gained 
by increasing the size of the elementary cell and studying the convergence of critical 
exponents. The computational effort then drastically increases and becomes compa- 
rable to that of the well known finite size scaling method based on the transfer matrix 
technique. The advantage of the RG method for hierarchical network models is that 
it yields qualitative aspects of the RG flow in phase space of the original model using 
only small scale numerics. 
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of the Deutsche Forschungsgemeinschaft . 
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